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Abstract. This article is about a problem in the numerical analysis of random operators. We 
study a version of the finite section method for the approximate solution of equations Ax — b in 
infinitely many variables, where A is a random Jacobi operator. In other words, we approximately 
solve infinite second order difference equations with stochastic coefficients by reducing the infinite 
volume case to the (large) finite volume case via a particular truncation technique. For most of the 
paper we consider non-selfadjoint operators A but we also comment on the self-adjoint case when 
simplifications occur. 

Mathematics subject classification (2000): 65J10; Secondary 47B36, 47B80. 
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1 Introduction 



Let U, V and W be non-empty compact subsets of the complex plane, and put 
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minM, w* := max lid and 5 := v* — (u* + w*). (I) 



We write N, Z, R and C for the sets of all positive integer, integer, real and complex numbers. 
Infinite matrices. In this paper we study bi- and semi-infinite matrices of the form 
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with entries Ui G U , Vi G V and w% G W for all i under consideration, where the box marks the 
matrix entry of A at (0, 0). As usual, we call u := (itj), v := (v{), and w := (wj) the sub-, main- 
and superdiagonal of A, resp. A + . We understand A and A + as linear operators, again denoted 
by A and A + , acting boundedly, by matrix- vector multiplication, on the standard spaces £ P (Z) 
and £ P (N) of bi- and semi-infinite complex sequences with p G [1, oo]. It is clear that the matrices 
(2) are in general not self-adjoint. We will study the selfadjoint case, when Wi = Ui+i for all i, 
separately in Section 2.6. 
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The sets of all operators A and A + from (2) with entries u, E U , v% E V and Wi £ W for all 
indices i that occur will be denoted by M(U,V,W) and M + (U,V,W), respectively. The set of 
all n x n matrices with subdiagonal entries in U, main diagonal entries in V and superdiagonal 
entries in W (and all other entries zero) will be called M n (U, V,W) for n E N, and we finally 
put M fin (C7, V, W) = U neN M n (U, V, W). 

Recall that a bounded linear operator B : X — > Y between Banach spaces is a Fredholm 
operator if the dimension, a, of its null-space is finite and the codimension, f3, of its image in Y 
is finite. In this case, the image of A is closed in Y and the integer ind^4 := a — (5 is called the 
index of A. For a bounded linear operator B on £ P (I) with I E {Z,N, Z \ N}, we write spec p B, 
speeds B and speCp t B for the sets of all A E C for which B — XI is, respectively, not invertible, 
not a Fredholm operator or not injective on £ P (T). Because A and A + in (2) are band matrices, 
their spectrum and essential spectrum do not depend on the underlying £ p -space [27, 29, 44], so 
that we will just write speci? and spec ess B for operators B in M(U, V, W) and in M + (U, V, W). 

Random alias pseudoergodic operators. Our particular interest is on random operators 
in M(U,V,W) and M + {U,V,W). We model randomness by the following concept: Given a 
metric space (Ai,d) and an index set I E {Z,N, Z \ N}, we say that a sequence a = (ctj)jgi i n 
A4 is pseudoergodic if for every e > 0, all n E N and all b = (bi)f =1 E M n , there is a k E I such 
that d{ak +i , b{) < e for all i = 1, n. In particular, if (A4 , d) is a discrete space then a = (aj) is 
pseudoergodic if and only if every finite vector over M. can be found (as a sequence of consecutive 
entries) in a. For a finite set M, a pseudoergodic sequence (ai)igN can be constructed by writing 
all -valued sequences of length 1, then 2, then 3, . . . in a row. For M = {0, 1}, this is done 
by stringing together the binary expansions of all natural numbers. 

We will call (ai)j e z : left-pseudoergodic if (a«)igz\N is pseudoergodic, right-pseudoergodic if 
(ctj)igN is pseudoergodic, and bi-pseudoergodic if both (ai)j e z\w and (ai)jgN are pseudoergodic. 
It is a simple exercise to show that (aj)jgz is pseudoergodic if and only if it is right- or left- 
pseudoergodic. 

Pseudoergodicity was introduced by Davies [13] to study spectral properties of random op- 
erators while eliminating probabilistic arguments. Indeed, if a = (aj)jgi, where all entries cij 
are independent (or at least not fully correlated) samples from a random variable with values 
(densely) in M then, with probability one, a is pseudoergodic in cases I E {N, Z \ N} and 
bi-pseudoergodic in case I = Z (e.g. [30, §5.5.3]). 

We call an operator A E M(U, V, W) pseudoergodic, left-pseudoergodic, right-pseudoergodic 
or bi-pseudoergodic and write A E VE(U,V,W), A E VE L (U,V,W), A E ^E R (U,V,W) or 
A E ^E 2 (U, V, W), respectively, if a = (af) i6 z with m := (ui, Vi, Wi) E M := U x V x W C C 3 
has the corresponding property. So we have 

9E(U,V,W) = VE L (U,V,W) U VE R (U,V,W), 
^E 2 (U,V,W) = VE L (U,V,W) n *Er(17,V,W). 

If A E 'I'Er^J/, V, W) then we will write A + E V I / E_|_(C/, V, W) for the corresponding semi-infinite 
submatrix A + of A from (2). We will say a little bit about spectral properties of pseudoergodic 
operators A and A + but will mainly focus on another problem: 

The finite section method (FSM). If one wants to solve an equation 

Ax = b, i.e. '^2 a ij x {j) = &0Oj ? E Z (3) 
jez 
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on X = £ p (7j) approximately, where A : X — > X (bounded) and b £ X are given and x G X 
is sought for, one often uses a projection method. Therefore, let P\ r : X — > X stand for the 
operator of multiplication by the characteristic function of the discrete interval Zn [/, ...,r] for 
/, r G Z with I < r and denote the image of Pi r by X\ r = C r ~ l+1 . One then picks sequences of 
integers I1J2, ■■■ — > —00 and r±, r%, ... — > +00 and replaces the infinite system (3) by the sequence 
of finite systems 

p in,r n Ap i n ,r n x n = Pi ntrn b, i.e. ^ Otf x n (j) = b(i), l n <i<r n (4) 

ln<j<r n 

with n G N . The aim is that, assuming invertibility of A (i.e. unique solvability of (3) for all 
b G X), also (4) shall be uniquely solvable for all sufficiently large n and the solutions x n G Xi n)Tn 
shall remain bounded in n and converge componentwise 1 to the solution x of (3). If the latter is 
the case for all right-hand sides del then we say that the finite section method (short: FSM) 
with cut-offs at (l n ) and (r n ) is applicable for A. 

If (ay) is a band matrix, i.e. dij = for \i — j\ > d with some d G N (as is the case for 
our operators (2)), then the FSM is applicable if and only if A is invertible and the sequence 
(A n ) := (Pi n ,r n APi n)rn ) n £H is stable [46]. By the latter we mean that there exists a no G N such 
that A n : Xi n>rn — > X\ n ^ n is invertible for all n > uq and sup n>no \\A~ \\ < 00. 

In the case l n = —n, r n = n we will speak of the full FSM for A. Recently, it has been shown 
in different situations that (and how) applicability of the FSM can be established by choosing 
the sequences (l n ) and (r n ) accordingly [32, 41, 42] if the full FSM is not applicable. We will 
formulate a condition for applicability of the FSM in the case of general sequences (l n ) and (r n ) 
and then apply this result to the case of pseudoergodic operators (2). For semi-infinite systems 
on X = £ P (N), replace Z by N and l n by 1 in all of the above. 

Motivation. A major motivation for the study of random Jacobi operators, their spectra 
and the solutions of the corresponding operator equations comes from condensed matter physics: 
Questions about the conductivity of certain (composed, disordered) media, about flux lines 
in superconductors or about systems of asymmetricly hopping particles have been modeled 
by random Schrodinger operators (Anderson model [1, 2]), non-self adjoint versions (Hatano & 
Nelson [21, 22, 23]) and other non-selfadjoint random Jacobi operators (Feinberg & Zee [16, 17]). 
Similar models arise in population biology [35]. Besides such discrete models also continuous 
problems that have been described by a stochastic differential equation in ID lead, after suitable 
discretization, to a matrix equation of the kind studied here. 

We give some upper and lower bounds on the spectrum of our operators but mainly focus 
on the approximate solution of operator equations Ax = b via the FSM. The latter can however 
be useful for spectral studies again: The inverse power method for the computation of the 
eigenvalue of A that is closest to a given point z G C approximates the (in modulus) largest 
eigenvalue of (^4 — zl)~ l by repeatedly solving equations (^4 — zl)x^ n+1 ^ = x^ n \ n = 0,1,..., 
with a rather arbitrary (non-zero) initial vector x^ . 

Historic remarks. The idea of the FSM is so natural that it is difficult to give a historical 
starting point. First rigorous treatments are from Baxter [3] and Gohberg &: Feldman [18] on 
Wiener-Hopf and convolution operators in dimension = 1 in the early 1960's. For convolution 
equations in higher dimensions N > 2, the FSM goes back to Kozak & Simonenko [25, 26], and 
for general band-dominated operators with scalar [37] and operator-valued [38, 39] coefficients, 

1 For p G (1, 00) this is equivalent [40] to convergence of the solutions x n (extended by zero) to x in X = £ P (Z). 
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most results are due to Rabinovich, Roch & Silbermann. For the state of the art in the scalar case 
for p = 2, see [45]. The quest for stable subsequences if the full FSM itself is instable is getting 
more attention recently [41, 42, 48, 49, 32]. In [42], the stability theorem for subsequences is 
used to simplify the criterion in dimension N = 1 by removing a uniform boundedness condition. 
However, we are not aware of a rigorous treatment of random (or pseudoergodic) operators via 
the finite section method. 

2 Main results 
2.1 Notations 

We first need some geometric notations: For sets S, T C C we put S + T := {s + t : s G S, t G T} 
and we write s + T := {s} + T and sT := {st : t G T} if s G C. By T = {z G C : \z\ = 1}, 
D = {z G C : \z\ < 1} and D = D U T we denote the unit circle, the unit disk and its closure. 
So, for example, S + eD is the closed e-neighborhood of S C C with e > 0. 




Figure 2.1: a) This is the ellipse E(u,w) with u = 3 and u; = i. The major axis of the ellipse bisects the angle 
between u and w at the origin. The half-axes (dotted lines) have length |3| ± |i|, i.e. 4 and 2. 
b) We see E+(U, W) in dark gray and E-(U, W) in light gray for U = {2} and W = {-1, 1}. E(U, W) is the 
union of the two ellipses E(2, -1) and E(2, 1). In Sections 2.2 and 2.3 we show that, for A £ *E([7, V, W) to 
be Fredholm (and hence invertible) it is necessary that either V C E+(U,W), in which case indA+ = —1, or 
V C E-(U, W), in which case xadA+ = 0. On the other hand, if V C D C £?+(?/, VF) then A is invertible. 

For u, w G C, put 

E(u,w) := {v€C : \v + 2y/^\ + \v - 2y/uw\ = 2(\u\ + \w\) } , (5) 

which is the ellipse that is centered at 0, has half-axes of length \u\ + \w\ and | \u\ — \w\ | and 
focal points ±2y/uw (so that the major axis of E(u,w) bisects the angle between u and w at 
the origin). By E + (u,w) and E-(u,w) we denote the bounded (interior) and the unbounded 
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(exterior) component of C \ E(u,w), respectively. Now, for non-empty U, W C C, let 



E(U,W) := [J E(u,w) and E±(U,W) := f] E±(u,w). 

ueu u eu 

iteff w g w 

Note that E(u, w) = -E(u,w) for all u,w G C, so that also E(U, W) = -E(U,W) and 
E±(U,W) = —E±(U,W) hold (see Figure 2.1 for an example). 



2.2 Spectrum and essential spectrum 

Let U, V, W C C be non-empty and compact sets, and recall (1). Then we have the following 
result about spectrum and essential spectrum of our pseudoergodic operators (2): 

Theorem 2.1 a) For A G VE{U, V, W) and A + G #E+(*7, V, W), the following holds: 
V + E(U,W) C 

|J spec B = (J spec^S = (J spec css B = [J spec css B+ 

BeM(uy,w) BeM(u,v,w) bgm(u,v,w) b + &m + {u,v,w) 

= spec ess A = spec A = spec ess A + C spec A + (6) 
C F+(n*+u;*)D 

b) The upper bound spec iC^| (u* + w*)D /rom (6) holds for arbitrary A G M(U, V, W) 
(as well as for semi-infinite and finite matrices A). For bi-infinite matrices A G M(U, V, W) we 
can improve this upper bound on spec A under one of the following conditions: 

if u* > w* , i.e. \u\ > \w\ Mu G U, w G W , then specj4 C C \ P| \v + (it* — w*)J^j, 

if u* < w*, i.e. \u\ < \w\ Mu G U, w G W , then spec AcC\ + (to* — ti*)B^ ; 

where, in addition to (1), we define := min^g^ |-u| and w* := min^gvv 1^1- 



We see from (6) that the spectrum of A and the essential spectrum of A and A + only depend 
on the sets U,V,W but not on the pseudoergodic operators A and A + . So all operators in 
| I / E(L r , V, W) have the same (essential) spectrum. In particular, in the case of random operators, 
spec ess A, spec A and spec ess A + do not depend on the distributions of the random variables for 
sub-, main- and superdiagonal - only on their supports U,V,W . (Note that none of the above 
applies to spec A + ; this set does depend on the concrete operator A + .) 

Example 2.2 a) Anderson model. In [1, 2], the conductivity of ID disordered media was 
studied. Here U = W = {1} and V C M, so that A is a discrete Schrodinger operator with 
random potential. In this case E(U, W) = [—2, 2]. So our lower and upper bound from Theorem 
2.1 are V + [—2, 2] and V + 2D. Together with spec A C R, by selfadjointness, we get that the 
lower bound is also an upper bound, whence all sets in (6) are equal to V + [—2,2]. 

b) Hatano 8z Nelson. The so-called non-selfadjoint Anderson model was introduced in 
[21, 22, 23] for the study of flux lines in type II superconductors under the influence of a tilted 
external magnetic field. Here U = {e 9 }, V = [—a, a] and W = {e~ 9 }, where a and g (the 
strength of the magnetic field) are positive real parameters. Now E := E(U, W) has half-axes 
of length e 9 ± e~ 9 being part of the real and imaginary axis, and E gets closer to a circle as 
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g — > oo. Abbreviate e 9 + e~ 9 = 2 cosh g =: c and e 9 — e~ 9 = 2sinh<7 =: s. If V is (at least) 
as long as the major axis of E, i.e. if a > c, then V + E and V + cB, which are the lower and 
upper bound in (6), only differ by 2e~ 9 in Hausdorff distance. It is easy to see that the closed 
numerical range of A (which is always an upper bound on spec A) is contained in V + conv(-E), 
which is equal to V + E if and only if a > c, so that all sets in (6) are equal to V + E in this case 
(cf. [12]). If a < c then the lower bound V + E has a hole around the origin, so that the best we 
can say then is V + E C spec A conv(£7). However, if V is shorter than the short axis of 

E, i.e. if a < s, then statement b) of the theorem proves that there is indeed a hole in spec A: 
Since = e 9 > e~ 9 = w*, we have that n veV (v + (it* - w*)B) = (-a + sB) n (a + sB) / 
is in the resolvent set of A. We summarize these bounds on spec^4 in case a < s in Figure 4.1 
below. A further study of the shape and size of the hole in spec A is in [12, 13, 14, 33, 34]. 

c) Feinberg & Zee. In [16, 17, 24] the case_C/ = {1}, V = {0}, W = T is studied. A 
simple computation shows that then E(U, W) = 2B = (it* + w*)B holds, so that all sets in (6) 
coincide. In the same papers the much more complicated case with W = {±1} is also studied. 
In this case, E(U, W) = [-2, 2] U [-2i, 2i] is far away from (u* + w*)D = 2B (see [6, 7, 8, 24] for 
sharper bounds in this case). □ 

Remark 2.3 The upper and lower bound in (6) might create the impression that spectrum 
and essential spectrum of A can be written as V + S(U, W) with a set S(U, W) independent of 
V, in which case it would be sufficient to study the case V = {0}. To see that this is not true, 
compare the cases UxVxW = {0} x_{±l} x {1} (see [50, 31]) and U x V x W = {0} x {0} x {1}: 
In the first case one has S(U, W) = B, whereas in the second case, S(U, W) = T. □ 

Both upper and lower bound in Theorem 2.1 can be improved: The lower bound comes from 
evaluating Uspec-B with the union taken over all B G M(U, V, W) that have constant diagonals. 
A better lower bound can be derived in concrete examples by also considering operators B € 
M(U, V, W) with diagonals of period 2, 3 or more (e.g. [6, 8, 12, 31]). The upper bound can be 
improved by different approaches such as (higher order) numerical ranges or hulls [14, 15] or by 
the more recent ideas of [7]. 

2.3 Fredholmness, invertibility and the full FSM 

Here are our results on invertibility, Fredholm property, and applicability of the full FSM. 

Theorem 2.4 If A G VE(U,V,W) or A + E WE+(U,V,W) is Fredholm then VnE(U,W) = 0. 
In fact, either 

(a) V^E^(U,W), or 

(b) V C E+(U, W) and u* > w* , or 

(c) V C E+(U, W) and u* < w*, 

where, in addition to (1), we define := mm U £u \u\ and w* := mm w ^\y \w\. 

The three cases correspond to the Fredholm index of A + (the so-called plus-index of A): 
(a) <^=> ind^ + = 0, (b) <==>- indA + = -1 and (c) <S=> mAA + = l. 

Note that, while the index of A + can be —1, or 1, the index of A is always zero if A is 
Fredholm; in fact, A is always invertible if Fredholm (see (6) or the following theorem). 
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Theorem 2.5 Let U, V, W C C be non-empty and compact. For A G M(U,V,W), we look at 
the following statements: 

(i) A is a Fredholm operator, 

(ii) A is invertible, 

(Hi) the full FSM is applicable to A, 

(iv) all operators in M (U, V, W), M + (U, V, W) and M fin (J7, V, W) are invertible 
and their inverses are uniformly bounded from above, 

(v) the full FSM is applicable to all operators in M(U, V, W), 

(vi) all B G M(U,V,W) are invertible, 

(vii) all B G M(U,V,W) are Fredholm operators. 

a) For general A G M(U, V, W), the following implications trivially hold: 

(i) <S= (ii) <S= (in) <S= 

1t it 1t (iv) 

(vii) 4= (vi) <= (v) ■<= 

b) If A £ fy~E(U,V,W) then (i), (ii), (vi) and (vii) are equivalent and (Hi) -<=^ (v) holds. 

c) If A G *$>E(U, V, W) and G U, W then (i) - (vii) are all equivalent. 

d) If 5 > in (1) then (i) — (vii) hold, where all the inverses are bounded above by 1/5. 

Remark 2.6 a) Statement b) of the theorem shows how 'hard' it is for a pseudoergodic 
operator to be Fredholm (i.e. invertible) or to even have an applicable full FSM. It also shows 
that, like the (essential) spectrum, these properties only depend on the sets U, V, W but not on 
the concrete operator A G ^E(U, V, W). 

b) If U, V, W are discrete sets and G U, W then A G ^E(U, V, W) decouples into a block 
diagonal operator A = Diag(£?j : i £ Z) with every B G M^ a (U, V, W) appearing as one of the 
blocks Bi, so that some of the above claims in c) become fairly obvious then. However, note 
that we do not assume U, V, W to be discrete in Theorem 2.5. 

c) The condition 5 > in d) is equivalent to V C C \ (u* + w*)U) or, to phrase it in the style 
of (6), to V + (u* + w*jB. □ 

Since applicability of the full FSM of an operator A G ^E(?7, V, W) is determined by the sets 
U, V, W only rather than by the operator, it seems advisable to use a version of the more flexible 
FSM (4) that gives credit to individual features of the concrete pseudoergodic operator A and 
will work under the sole condition of invertibility of A, where the full FSM might fail. 

We say 'might fail' because we do not have an example of sets U, V, W and A G ^/E(U, V, W), 
where the full FSM fails while the following version applies, unless when indA + ^ 0. However, 
we can prove that our adapted FSM from Section 2.4 generally applies if A is invertible - which 
we doubt in the case of the full FSM (even if ind^4+ = 0). 

2.4 The FSM with adaptive cut-off intervals 

If Theorem 2.5 does not yield applicability of the full FSM for A G ^>E(U, V, W) then we propose 
using the FSM (4) with cut-offs at integer values (/„) and (r n ) that are adapted to the operator 
A at hand. 



7 



The adaptive FSM in the general case. We start with a statement for general tridiag- 
onal operators A G M(U, V, W) or, in fact, for even more general operators. For X = £ P (Z) with 
p G [l,oo], we write A G BO(X) and call A a band operator if ^4 acts via matrix- vector mul- 
tiplication by a band matrix. Moreover, we write A G BDO(X) and call A a band-dominated 
operator if A is the limit (in the operator norm induced by || • \\x) of a sequence of band operators. 

If A G BDO(X) is given by the matrix (aj,-)i,-gz and -B G BDO(X) is given by a matrix 
(6ij)i,jez then we call i? a Jimit operator of ^4 if there exists a sequence h = (hi, hi, ...) of integers 
with \h n \ —> 00 and ai+h n ,j+h„ —> hj as n — > 00 for all z,j G Z. In this case we write .B =: ^4^. 
For a given sequence of integers going to infinity, let 

a k n3 (A) := {A g : 5 is an infinite subsequence of h for which ^4 5 exists}. 

By a Bolzano- Weierstrass argument it can be seen that a k v (A) is always nonempty. We will also 
abbreviate 

<j op (A) := } (A) and a°_ p (A) := a^^ JA) 

and put a op (A) := a°^(A) U a°^(A), so that the latter is the set of all limit operators of A. 

In the semi-infinite case X = £ P {H), the spaces BO(X) and BDO(X) are defined in the same 
way. For A G BDO(X) one then also defines limit operators Ah exactly as above - provided 
that h = (hi, hi, ...) tends to +00. Note that Ah is in any case bi-infinite, i.e. it acts on l v (7L). 

The set of pseudoergodic operators can be equivalently characterized in terms of limit oper- 
ators (see [28, §3.4.10] or [30, §5.5.3]): 

Lemma 2.7 For an operator A G M(U,V,W), one has 

( VE h (U,V,W) <=> a°_ p (A) ) 
A G I *E R (U,V,W) ^ a° + p (A) \ = M(U,V,W). 
[ ^E(U,V,W) a op (A) ) 

Let X = £ P (Z) and P : X — > X denote the operator of multiplication by the characteristic 
function of N, and let Q := I — P be the complementary projector of P. Given an operator 
A G BDO(X) with matrix (dij)ijez, we write A + for the compression PAP\ ira p of A to imP = 
^ P (N); that is, A+ is the operator of multiplication by the matrix (aij)ijeN- Analogously, we 
write A- for the compression QAQ\{ m Q of A to imQ = £ P (Z \ N); that is, A_ is the operator of 
multiplication by the matrix (tJy)i,j'eZ\N- When talking about their invertibility, Fredholmness 
or index, we always understand A + and A- as operators on £ P (N), resp. £ P (Z \ N). If we are 
only interested in an operator on £ P (N), we usually denote it by A + (indicating that it is the 
compression of an operator A on -P(Z) to the positive half-axis) to remind ourselves of the 
semi-infinite setting. 

The following theorem is a generalization of results from [42, 32] (which can be derived by 
straightforward changes in the proofs there): 

Theorem 2.8 Let X = £ P (Z) with p G [1, 00] and fix two sequences I = (Z n )raeN and r = ( r n)neN 
of integers I1J2, ■■■ — > —00 and r\,ri, ■■■ —> +00. For A G BDO(X), the finite section method 
(4) is applicable if and only if the following operators are invertible: 

A, all operators B + with B G o~° p (A), all operators C_ with C G a° p (A). (7) 
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The set of operators (7) is particularly handy if the sequences I and r are such that 
the sets {B + : B £ cr° p (A)} and {C_ : C £ cr° p (A)} are singletons, 
which is equivalent to the existence of the strong limits B + and C_ of 

and 



(8) 



U l n +1 vi n + i w ln + 1 

Uln+2 V in+2 



\ 



\ 



Vr n -2 Wr n -2 
u r n -l «m-l w r n -l 



as n — > oo. For (8) it is sufficient (but not necessary) that the limit operators A[ =: B and 
A r =: C exist. 

Here is the version of Theorem 2.8 for semi-infinite matrices: 

Theorem 2.9 Let X = £ P (N) with p G [l,oo] and fix a monotonously increasing sequence 
r — ( r n)neN of positive integers. For A + £ BDO(X), the finite section method (4), toitft Z n = 1 
/or aZZ n G N, is applicable if and only if the following operators are invertible: 



A + , all operators C_ with C £ a° p (A + ). 
Also here, the set (9) is smallest possible if {C_ : C G o~° p (^4+)} is a singleton. 



(9) 



Bi-infinite bi-pseudoergodic systems. We demonstrate how, under the sole (and for this 
purpose minimal - because necessary) assumption of invertibility of A, one can approximately 
solve operator equations Ax = b on £ P (Z) with a bi-pseudoergodic operator A by the finite 
section method. 



Algorithm 2.10 — The ^ / E2-FSM. Suppose U,V,W C C are non-empty and compact sets, 
p G [1, oo], A G *E 2 (C7, V, VP) is invertible and b G F(Z) is given. 

Step 1. Pic& some arbitrary u £ U , v £ V and w £ W. Choose integer sequences h,l2,~- 
monotonically decreasing and r%,r2, ■■■ monotonically increasing such that l n < r n and 



1 

\Ui — U\ + \Vi — V\ + \Wi — w\ < — , 

n 



Vi G {l n , l n + 1) — , l n +n}U { v n n, . . . , r n 



forn = 1,2,... . 

Step 2. By Theorem 2.4, we know that we are in one of the three cases (a), (b), (c). To 
find out which of these cases applies, compute 



> 



\v + 2\Juw\ + \v — 2\/uw\ — 2(\u\ + 



case (a), 



< : compute \u\ — \ w\ 



> 
< 



case (b), 
case (c). 



If one of the expressions to be computed here is zero or if the outcome of this algorithm depends 
on the choice ofu,v,w then A is not Fredholm, let alone invertible, by Theorem 2.4- 

Step 3. Depending on case (a), (b) or (c), we apply our finite section method (with cut-offs 
at (l n ) and (r n ) as chosen in step 1) to different equations: 
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In case (a), the FSM (4) is applicable to the equation Ax = b, i.e. to 



t>_ 2 W-2 

U-l V-l W-l 

u v w 



Ul 



V\ Wl 
U 2 V 2 



x(-2) 
x(-l) 
x(0) 
x(l) 
x{2) 



6(-2) 
fe(-l) 

Kg) 

6(1) 
6(2) 



In case (b), the FSM (4) is applicable to the following upward-translated (obviously equivalent) 
system 

(■.-.-. \ ( : \ / : \ 



V-i W-l 



Uq Vq w 



U\ V\ Wl 
U2 V2 

a?, 



V 



x(-2) 
x(-l) 
x(0) 
x(l) 
x(2) 



V 



K-i) 
&(o) 

6(1) 
6(2) 
6(3) 



V 



Finally, in case (c), the FSM (4) is applicable to the following downward-translated (obviously 
equivalent) system 



( : 



\ ( 



W-3 

U_l W-l 



"0 



Vq Wq 



Vl Wl 



■ / 



x(-2) 
x(-l) 
x(0) 
x(l) 
x(2) 



I 



\ 



6(-3) 
6(-2) 

&(-!) 
6(0) 
6(1) 



Remark 2.11 — The growth of the intervals {l n , r n } as n — >• oo. 

a) If all entries Ui, V{ and ujj of A are independent samples from three uniformly distributed 
random variables with values (everywhere) in U, V and W then, for all choices u, v, w in step 
1, one expects the same exponential growth of —l n and r n . For example, if U,V,W are finite 
with \U x V x W\ = m then —l n and r n are of order m n . (In the random but not uniformly 
distributed case, it is certainly advisable to pick some of the more likely u, v, w in step 1 in order 
to minimize the growth of —l n and r n .) 

b) The choice of l n and r n in step 1 is such that the conditions of Theorem 2.8 and in 
particular condition (8) are met with operators B + and C_ having constant diagonals (containing 
u,v and w). It is possible to aim at different (non-Toeplitz) operators B + and C_ here (as long 
as they are invertible) via the choice of l n and r n , while, possibly, keeping the growth of —l n 
and r n more moderate. Steps 2 (with arbitrarily picked u, v, w) and 3 still remain as shown. 

c) One should not be too worried if these finite systems become large very quickly since 
they can be solved in linear time (as opposed to cubic, for the Gauss algorithm). In case (a) 
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this is done by the so-called Thomas algorithm [11], while in case (b), resp. (c), the solution is 
calculated successively via backward, resp. forward, substitution. 

d) In [41] the finite section method (4) is adapted to the Almost Mathieu operator 

{Ax) n = x n _i + A cos (2ir(na + 8)) x n + x n+ i, n£Z 

by putting — l n = r n equal to the denominator of the n-th continued fraction approximant of the 
irrational number a G (0, 1). Note that this means —l n = r n also grow exponentially in n. For 
example, if a = (v5 — 1) /2 is the golden mean then —l n = r n is the n-th Fibonacci number. The 
order of exponential growth is higher if the continued fraction expansion of a contains larger 
numbers. (For the golden mean, it is 1/(1 + 1/(1 + 1/ • • • )).) □ 

Semi-infinite pseudoergodic systems. For semi-infinite systems A + x = 6 on £P(N), the 

situation is related but much simpler. Again, we only assume invertibility of the operator. 

Algorithm 2.12 — The ^E + -FSM. Suppose U, V, W C C are non-empty and compact sets, 
p G [l,oo], A + G tyE + (U,V,W) is invertible and b G £ P (N) is given. 

Step 1. Pick some arbitrary u G U , v G V and w G W . Choose a monotonically increasing 
sequence ri,r2,... of positive integers such that 

1 

\Ui — U\ + \V i — V\ + \Wi — w\ < — , 



Vi G {r n - n, ...,r n - l,r„} 
it 

forn = 1,2,... . 

Step 2. By Theorem 2.4 and the invertibility of A + , we are automatically in case (a). 

Step 3. The FSM (4) with l n = 1 for all n G N and (r n ) as chosen in step 1 applies to our 
equation A + x = b, i.e. to 

( v x wx \ ( x(l) \ ( 6(1) \ 



U 2 V2 W2 

u 3 v 3 w 3 

U4 «4 



V 



/ 



x{2) 
x(3) 
x(4) 

V = J 



6(2) 
6(3) 
6(4) 

V '■■ J 



As in Remark 2.11 b), note that one could choose r = (r n ) in step 1 so that C_, with 
C G <7° P (A), is not of Toeplitz structure but is another (invertible) operator. For example, one 
could choose r„ such that 



\Ur n -i-Ui + 2\ + \Vr n -i-V i+1 \ + \W rn - 
for n = 1,2, so that 



Wj\ < 



1 



n 



Mi G {0, ...,n} 



/ 



V 



-2 W rn - 2 

-1 Wfn-1 W r n 



\ 



J 



\ 



\ 

v 3 w 2 

U 3 V 2 Wl 

U 2 Vi J 



C- 



(10) 



strongly as n — > oo. But C_ is invertible by our assumption on A + : We call the matrix C_ in 
(10) the reflection of the operator A + from (2) and we will write A^ for C_. Conversely, we also 
call A_|_ the reflection of C_ and denote it by C^. It is easy to see that a semi- infinite matrix is 
invertible if and only if its reflection is invertible. 
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2.5 Spectral and pseudo-spectral approximation 

Here we will briefly discuss another feature of the finite section method with adaptive cut-off 
intervals. We will work exclusively on the Hilbert space H := £ 2 (Z,). Let again I and r be 
sequences of negative and positive integers which converge to — oo and +00, respectively. The 
set Ai >r of all band-dominated operators A on H for which the limit operators A\ and A r exist 
is a C*-algebra, as one easily checks. Note that every band-dominated operator A belongs to 
an algebra Ai, r with specified sequences l,r. 

Proposition 2.13 Let A G Ai >r . Then the adaptive finite section method (4) with system 
matrices Pi n ,r n APi n>rn is fractal. 

The notion of a fractal approximation method was introduced in [47]. Since already its defi- 
nition makes heavily use of C*-algebraic language we will omit all technical details here and 
refer the interested reader to [47] and [20]. Roughly speaking, an algebra of approximation 
sequences is fractal if every sequence in the algebra can be reconstructed from each if its (in- 
finite) subsequences modulo a sequence which tends to zero in the norm. A single sequence 
like (Pi n ,r n APi n) r n ) n ^ is called fractal if the smallest C*-algebra which contains this sequence 
and the sequence (P/ njr . n ) ng pj has the fractal property. The proof of Proposition 2.13 follows 
easily from Theorem 2.8 above and Theorem 1.69 in [20]. The main point is that the sequence 
(Pin,r n APi n ,r n )n&i is stable by Theorem 2.8 if and only if the operators A, B + = PAiP and 
C- = QA r Q are invertible and that the operators A, B + and C- are strong limits of (shifts of) 
the sequence {Pi n ^ n APi n:rri ) n< ^. Since every subsequence has the same strong limits, the result 
follows from Theorem 1.69 in [20]. 

Fractal sequences are distinguished by their excellent convergence properties. To mention 
only a few of them, let a(A) denote the spectrum of an operator A, write 02 (^4) for the set of 
the singular values of A, i.e., o~ 2 (A) is the set of all non- negative square roots of elements in the 
spectrum of A* A and finally, for e > 0, let o~( £ \A) refer to the e-pseudospectrum of A, i.e. to 
the set of all A G C for which A — XI is not invertible or \\(A — > 1/e. Let further 

dn(M, N) := maxjmax min \m — n\, max min \m — n\} 

meAf neW n£N m<=M 

denote the Hausdorff distance between the non-empty compact subsets M and N of the complex 
plane. 

Theorem 2.14 Let A G A^ r and A n := P\ n ,r n APl n ,r n - Then the following sequences converge 
with respect to the Hausdorff distance as n — > 00 : 

(a) a{A n ) a(A) U a{B + ) U a(C_) if A is self-adjoint; 

(b) a 2 (A n ) -^a 2 (A)Ua 2 (B + )Ua 2 (C.); 

(c) a^{A n ) o-^iA) U a^(B + ) U a^(C^). 

The proof follows immediately from the stability criterion in Theorem 2.8, from the fractality of 
the sequence (A n ) by Proposition 2.13, and from Theorems 3.20, 3.23 and 3.33 in [20]. Let us 
emphasize that in general one cannot remove the assumption A = A* in assertion (a), whereas 
(c) holds without any assumption. This observation is only one reason for the present increasing 
interest in pseudospectra. For detailed presentations of pseudospectra and their applications as 
well as of other spectral quantities see the monographs [4, 5, 20, 51] and the references therein. 
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2.6 The selfadjoint case 

We discuss briefly how our results simplify when A is selfadjoint, i.e. w% = Th+i for all % in (2). 

In that case, there are only two sets, U and V, of which V is real. Instead of the ellipses 
E(u,w), one looks at E(u,u) = [— 2\u\, 2\u\]. The set E(U,W) gets replaced by the union of 
E(u,u) over all u 6 U, which is simply [— 2u*,2u*]. 

In Theorem 2.1 a), the lower bound therefore becomes V + [—2u*, 2u*\. But the upper bound 
becomes the same (recall Example 2.2 a) since w* = u* and since all spectra are real. So 

spec^4 = spec css vl = spec ess ^4 + = spec^4 + = V + [— 2u*, 2u*]. (11) 

In particular, A is positive definite if and only if min^y > 2u*. 

In Theorem 2.4 and anywhere else, case (a) applies (all indices are zero of course). The 
theorem says that V and [— 2u*,2u*] are disjoint if A is Fredholm. From (11) we know that 
V n [— 2u*, 2u*\ = is indeed both necessary and sufficient for A to be Fredholm (i.e. invertible). 

Concerning the FSM, not much simplification occurs apart from the fact that the full FSM 
is applicable if A is positive or negative definite. 

3 Background theory and proofs 

Before we come to the deeper results, let us briefly show how Fredholmness (and index) of A 
is related to that of its half-axis compressions A+ and A-. The following lemma is taken from 
[36, 43]. 

Lemma 3.1 An operator A £ BO{X) is Fredholm if and only if both its compressions A + and 
A- are Fredholm, i.e. spec ess A = spec ess A + U spec ess A- . Moreover, ind A = ind A + + ind A- . 

Proof. Since PAQ and QAP are of finite rank if A is a band operator, one has that A = 
PAP + PAQ + QAP + QAQ is equivalent, modulo compact operators, to 

PAP + QAQ = {PAP + Q)(P + QAQ) = (P + QAQ)(PAP + Q). 

So A is Fredholm if and only if PAP + Q and P + QAQ, which are the extensions (by identity) 
of A + and A- to -P(Z), are Fredholm. Moreover, 

ind A = ind (PAP + PAQ + QAP + QAQ) = ind (PAP + QAQ) 

= \nd(PAP + Q)(QAQ + P) = ind (PAP + Q) + ind (QAQ + P) 
= ind A + + ind A_ 

holds. ■ 

One refers to indj4_ and indj4 + as the minus- and the plus-index of A. By Lemma 3.1, the 
problem of determining Fredholmness (and the index) of A splits into two subproblems. These 
two subproblems again split into many smaller problems, where the key notion is again that of 
a limit operator. Besides Lemma 2.7 and Theorems 2.8 and 2.9, limit operators feature in the 
following characterization of Fredholmness (including the index) : 
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Theorem 3.2 Let X = £P(I) with p £ [1, oo] and I € {Z, N, Z \ N}, and let A £ BO(X). 

a) T/ie following are equivalent 

(i) A is Fredholm on X, 

(ii) all limit operators of A are invertible on £ P (Z) [37, 40], 
(Hi) all limit operators of A are injective on £°°(Z) [9, 10], 

so that, by applying the above to A — XI in place of A, 

spec ess ^ = (J specS = |J spec^S. (12) 

B£a°P(A) Becr°P(A) 

b) If A £ BO(£ p (7*)) is Fredholm then all operators in o~°^(A) have the same minus-index 
and all operators in a°^(A) have the same plus-index, which also happen to be the minus- and 
the plus-index of A, respectively [36, 43]- This means ind A = indj4_ + ind^4 + , where 

indA_ = mdB^ V B G a op (A), (13) 
and mdA + = indC+ V C £ a°^(A). (14) 

Let S denote the shift operator (Sx)(m) := x(m— 1), m £ Z, on X = £ P (Z). If A is pseudoergodic 
then a op (A) = M(U,V,W), by Lemma 2.7. Particularly simple elements of M (U, V, W) are 
operators whose matrix has constant diagonals. So fix u £ U, v £ V and w £ W, and let 
L(u, v, w) := uS + vl + wS" 1 be the single element of M({u}, {v}, {w}) C M(U, V, W), which 
is a so-called Laurent operator (sometimes also called "bi-infinite Toeplitz operator"). It is a 
standard result [4, 5] that 

spec L(u,v,w) = {ut 1 + vt° + wt' 1 : t £ T} = v + E(u,w). (15) 

Together with (12), the latter proves the lower bound in Theorem 2.1. The rather crude (but 
still helpful) upper bound in Theorem 2.1 and statement d) in Theorem 2.5 rely on the following 
simple lemma and its corollary. 

Lemma 3.3 Let A be a finite or (semi- or bi-)infinite matrix with subdiagonal u = (uj), main 
diagonal v = (vi) and superdiagonal w = (wi). Put 



u* : = sup|uj|, := inf|uj|, w* := sup|u;i|, and 5a '■= v* — (u*+w*). (16) 

i 1 i 

If & A > then A is invertible and H^" 1 ]] < 1/6 a- 

Proof. Write A = D + T with D = diag(uj) and treat A as a perturbation of D. We have 
A = D{I + D _1 T), where D _1 T has subdiagonal entries lijf" 1 , superdiagonal entries WiV^ 1 and 
everything else zero. From 5a > we get that 

HT^-lmii -li -li u* + w* 

\\D T\\ < supp^ | + sup\WiV i | < < 1, 

i i ^* 

so that I + D~ 1 T is invertible by Neumann series. But from A" 1 = (I + D^T)- 1 D' 1 it follows 
that also A is invertible and 

iia-ii < iKz+K-mniD-'ii < I - p L_^| C -i|| < TI 4 Sf i = 

II II V* ^ 

as was claimed. ■ 
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Corollary 3.4 IfU,V,W C C are non-empty and compact, (1) holds with 5 > 0, and if A is in 
M(U,V,W) orM + (U,V,W) or M^{U,V,W) then A is invertible and < 1/5. 

Proof. Just note that 5a > 5 for 5a from (16) and 5 from (1) and apply Lemma 3.3. ■ 

Now we have all the machinery to prove our main results: 

Proof of Theorem 2.1. a) All unions in this proof are taken over the set of all B G 

M(U,V,W). By (12) and Lemma 2.7, 

spec ess A = U spec^ B = U spec B D U spec ess -B 5 spec ess ^4, 

so that equality holds in both "5" signs. Moreover, 

spec ess A C spec j4 C U spec £> = spec ess ^4 

holds since A is one of the operators B in this union. So again we have equality everywhere. 
Equality (12) also holds with A replaced by A+. Hence, by cr op (A + ) = M(U,V,W), 

spec A D spec css A D s pec ess A + = U spec B D spec ^4 

holds, which proves the remaining equality in (6). The lower bound V + E(U,W) in (6) now 
follows by evaluating specB from (15) for all Laurent operators B = L(u,v,w) G M(U,V,W). 
The upper bound V + (it* + w*)D follows from Corollary 3.4 since A + — XI + is invertible if 
\v — A| > |u| + \w\ for all (u, v, w) G U x V x W, i.e. if dist(A, V) > u* + tM*. 

b) Let A £ M(U, V, W), u* > w* and suppose A G (^^£^(1; + («* - w*)0). Then |« - A| < 
u* — w* < \u\ — \w\ for all u G J7, u G V and it; G W, so that the subdiagonal of j4 — A/ 
dominates the other two diagonals. By a simple perturbation argument as above (see Lemma 
3.3 and Corollary 3.4), S'~ 1 (A — XI) = diag(itj)(7 + T) with ||T|| < 1 is invertible, and hence 
A — XI is invertible. The argument for the case w* > u* is completely symmetric. ■ 

Proof of Theorem 2.4. If A is Fredholm then all its limit operators B, including the Laurent 
operators B := L(u,v,w) G M(U,V,W), are invertible. So, for all (it, v, w) G U x V x 14 7 , we 
have that specL(it, u, id) = u + E(u, w) = v — E(u, w), i.e. v G" E(u, w). The following three 
cases are possible: 

(a) wind(-E(u, w), v) = 0, i.e. v is in the exterior of the ellipse E(u,w), or 

(b) wmd(E(u,w),v) = 1, i.e. v is encircled counter-clockwise by E(u,w), or 

(c) wind(E (u,w),v) = —1, i.e. v is encircled clockwise by E(u,w), 

where wind(C, z) denotes the winding number of a closed oriented curve C w.r.t. a point 
z G 1 C and where the ellipse E(u, w) is parametrized (and thereby oriented) by the map <p 1— >■ 
ite 1 ^ + we~ l{p from [0, 27r) to E(u,w). A simple computation shows that E(u,w) is oriented 
counter-clockwise if |it| > \w\ and clockwise if |u| < \w\ (while the ellipse degenerates into a line 
segment if |u| = \ w\). Let g denote the rotation z *— > | — z of the complex plane around |. For 
the Toeplitz operator B + , one has (see e.g. [4, 5]) 

indi? + = — wind(spec B, 0) = — wind(i> + E(u, w), 0) = —wmd(g(v + E(u,w)),g(0)) 
= — wind(— I — E(u, w), |) = — wind(— E(u, w), v) = —w'md(E(u,w),v), 

which is in case (a), —1 in case (b) and 1 in case (c). By (14), we have indi?+ = ind^4 + for 
all B G ct° p (j4) = M(U, V, W), so that for all choices (it, v, w) G U x V x W, the same case, (a), 
(b) or (c), applies - according to ind^4+. ■ 
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Proof of Theorem 2.5. a) The only implication that is not obvious here is that A is invertible 
if the full FSM of A is stable. This can be found in [46] (also see [40, 28]). 

b) Now let A G ^E(U,V,W). Properties (i), (ii), (vi) and (vii) are equivalent because (i) 
implies (vii) by Theorem 3.2 a) and Lemma 2.7. It remains to show that (Hi) implies (v). By 
Theorem 2.8 with I = (—1, —2, ...) and r = (1, 2, ...), property (Hi) is equivalent to invertibility 
of A and all operators B + and C_ with B G a 09 (A) and C G a°^(A). W.l.o.g suppose A is 
right-pseudoergodic, so that a°^(A) = M(U, V,W) by Lemma 2.7 and hence C_ is invertible 
for all C G M(U, V, W). But then, for every B G M(U, V, W), also B + is invertible because its 
reflection B+ is of the form C_ for some C G M(U, V, W) and is therefore invertible. Finally, 
since every operator in M(U,V,W) is invertible if A is invertible (see above), we conclude (v), 
by Theorem 2.8 again. 

c) Now let A G ^E(U, V, W) and G U, W. To see that all properties (i)-(vii) are equivalent, 
it is sufficient to show that (i) implies (iv). So let A be Fredholm. By Theorem 3.2 a) we 
know that all limit operators B of A, which are all operators in M(U,V,W) by Lemma 2.7, 
are invertible. Moreover, there is a c > (e.g. the norm of a Fredholm regularizer of A, 
see [37]) such that ||-B _1 || < c for all these operators B. Now we can show (iv): If D G 
M(U,V,W) then B := D G M(U,V,W) is invertible and WD" 1 ]] < c. If D + G M + (U,V,W) 
then B := Diag(L>^, L>+) is in M(U, V, W) since G U, W and hence B is invertible, so that D + 
is invertible and ||( J D + )- 1 || = < c. Finally, if D G M fin (J7, V, W) then, since G U,W, 
B := Diag(--- ,D,D,D,---) G M(U,V,W) is invertible, so that D is invertible and ||£> _1 || = 
WB- 1 ]] < c. 

d) If 5 > then, by Corollary 3.4, property (iv) holds and hence all the others follow. ■ 

Proof of correctness of Algorithm 2.10. The choice of the sequences I = (l n ) and r = (r n ) 
in step 1 is such that the sets {D + : D G (J° p (^)} =: {B + } and {£>„ : D G cr° p (yl)} =: {C_} 
are singletons; in fact, B + and C_ are Toeplitz operators with diagonals u, v, w. It is possible 
to choose sequences I and r with these properties because A G ^^(U, V, W). 

The test in step 2 exactly follows the geometric definition (5) of the ellipse E(u,w): If 
\v + 2\Juw\ + \ v — 2y/uw\ > 2(\u\ + \ w\) then v G E_(u, w) and we are in case (a) of Theorem 2.4 
(also see the proof of Theorem 2.4). If \v + 2y/uw\ + \v — 2y/uw\ < 2(\u\ + \w\) then v G E + (u, w) 
and it remains to check the orientation of the ellipse. For \u\ > \w\, the ellipse is counter- 
clockwise oriented, so that case (b) applies, and for |ii| < the orientation is clockwise and we 
are in case (c). By Theorem 2.4, the outcome of this test does not depend on the values u G U, 
v G V, w G W chosen in step 1 if A is Fredholm. The resulting case corresponds to ind^4+. 

If we are in case (a) then indA + = 0. Because A is invertible, we have = indA = 
ind^4 + + ind^4_ = ind A—. Now let D G a° p (A). By Theorem 3.2 a) and b), D is invertible 
and ind-D+ = indA + = 0, so that = ind-D = indL> + + ind-D_ = indj4+ + ind(7_ = indC_. 
So C- is a Toeplitz operator that is Fredholm with index 0. By Coburn's theorem [4, 5], C_ is 
invertible. By a completely symmetric argument (or simply by noting that B + = C3) we get 
that also B + is invertible. Since (8) holds, Theorem 2.8 yields the applicability of the FSM (4) 
with the sequences / and r as chosen in step 1. 

If we are in case (b) or (c) then k := ind^4 + = =f1 and there are no sequences I and r for 
which the FSM (4) could be applicable to Ax = b. However, the FSM (4) is applicable to the 
equivalent system S k Ax = S k b by the same arguments as in case (a) since S k A is invertible and 
since m.&(S k A) + = ind S%A + = indS^ + ind ,4+ = -k + k = 0. ■ 
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Proof of correctness of Algorithm 2.12. The choice of the sequence r = (r n ) in step 1 is 
such that the set {D_ : D £ a° p (A)} =: {C_} is a singleton; in fact, C_ is a Toeplitz operator 
with diagonals u, v, w. It is possible to choose such a sequence r because A+ € \I/E_|_([/, V, W). 

Since ind ^4+ = by invertibility of A+, we are automatically in case (a) of Theorem 2.4. 
Let D 6 a° p (A). By Theorem 3.2 a) and b), D is invertible and indD + = mdA + = 0, so that 
= ind D = ind D + + ind Z)_ = ind A + + ind C_ = ind C_ . So C_ is a Toeplitz operator that is 
Fredholm with index 0. By Coburn's theorem [4, 5], C_ is invertible. Now Theorem 2.9 yields 
the applicability of the FSM (4) with the sequence r as chosen in step 1. ■ 

4 A numerical example 

We illustrate our results by a numerical computation, for which we come back to the Hatano- 
Nelson model from Example 2.2 b). So let U = {e 9 } and W = {e~ 9 } with g > 0, put c := 
e 9 + e~ 9 = 2 cosh g and s := e 9 — e~ 9 = 2 sinh g, and let V = [—a, a] with < a < s < c. 

Now let A £ ^~E(U, V, W). From Theorem 2.1 and our discussion in Example 2.2 b), we 
derive the upper and lower bounds on spec A as shown in Figure 4.1. For further studies of this 
operator, including the size and shape of the hole in its spectrum, see [12, 13, 14, 33, 34]. 




-c-a -s-a a-c a-s s-a c-a s+a c+a 

-c-a a c 



Figure 4.1: Here are our lower (dark gray) and upper (dark+light gray) bound on spec A from Example 2.2 b) 
with g = 1 and a — 2, so that < a < s < c. The region of uncertainty (light gray) is small if s and c are close to 
each other (i.e. if g is large). For reasons of symmetry we have only shown the upper half of the complex plane. 

While the spectrum of a n-by-n principal submatrix of A is less interesting (each such matrix 
is similar to a self-adjoint matrix), the spectrum of this finite problem with periodic boundary 
conditions and its limit as n — > oo has been described in much detail by Goldsheid and Kho- 
ruzhenko [19], who thereby verified numerical observations of Hatano and Nelson [21, 22, 23]. 
The limiting set as n — > oo turns out to be the union of certain analytic curves (the so-called 
'bubble with wings' [50]) and is entirely different from (although contained in) the spectrum of 
the infinite matrix A. 

We will now apply our adaptive FSM (Algorithm 2.10) to a concrete matrix A £ M(U, V, W), 
whose main diagonal entries V{ have been chosen independently from V = [—a, a], where the 
density of our probability distribution on V increases in a certain way towards the endpoints 
of the interval. Our model has the parameters g = 1 (so that c = 2 cosh 1 ~ 3.0862 and 
s = 2 sinh 1 ~ 2.3504) and a = 2 < s, whence spec A. 

In step 1 of the algorithm, we choose, as motivated in Remark 2.11 a), v = 2, which is one of 
the values with the highest probability density, besides the obvious choices u = e 1 and w = e _1 . 
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Then, for n = 1, 2, we look for n consecutive entries of the main diagonal that are within 1/n 
of v = 2 to find our cut-off bounds l n < l n -i and r n > r n _i (see Figure 4.2) with Iq = = rrj. 




-50 -40 I -30 -20 1-10 I r r 10 20 30 r 40 50 

3 2 1 1 2 3 

Figure 4.2: These are the main diagonal entries w_so to vso that are close to 2. Encircled are the groups of 
n = 1, 2, 3 consecutive entries that are within 1/n of « = 2 and therefore lead to the definition of l n and r n . 



In step 2 of the algorithm, we find that we are in case (b), which says that v lies inside 
the ellipse E = E(u, w) and it is encircled counter-clockwise w.r.t. the parametrization (p \— > 
ue iip + we - '*' = e 1+iip + e -1- ^ of E. In other words: ind^+ = -1. 

This means that, in step 3, we shift our infinite system up by one row before we truncate 
it according to our sequences I = (l n ) and r = (r n ). The resulting method is applicable if and 
only if the inverses of the finite matrices A n := Pi n ,r n S~ l APi n>rn remain uniformly bounded as 
n — > oo. The following table shows the cut-off sequences / = (l n ) and r = (r n ), the size of the 
matrices A n and the norms of their inverses for n = l,2,...,8. 



n 




r n 


f"n In ~\~ 1 


II A~ l \\ 
\\ A n II 


1 


-1 


1 


3 


0.6816 


2 


-12 


5 


18 


1.0580 


3 


-35 


35 


71 


1.2698 


4 


-41 


162 


204 


1.2698 


5 


-899 


537 


1437 


1.4121 


6 


-1068 


1183 


2252 


1.5438 


7 


-20494 


21758 


42253 


1.6135 


8 


-469241 


41570 


510811 


1.7500 



We see the rather irregular exponential growth of the intervals {l n ,...,r n } (see Remark 2.11) 
and the moderate growth of the inverses A^ 1 . This numerical evidence is not really convincing 
that the inverses remain uniformly bounded as n — > oo. However, from the theory behind our 
Theorem 2.8 it follows that lim sup n H^.^ 1 1| is in case p = 2 equal (and otherwise at least bounded 
above by two times) the maximum of the norms of the inverses of the operators in (7). In our 
case, this means that 

limsup ||^4^|| = max( \\A _1 \\ , 1| ) 

n 

if p = 2, where B + is the Toeplitz operator (note the translation S~ l in step 3) 

/ e 1 2 e" 1 \ 
e 1 2 ■•■ 

B + = 

e 

V '• / 

with symbol a(t) = e 1 +2t' 1 +e~ 1 t~ 2 , t G T. But from < (s-2) -1 « 2.8539 (the argument 

is as in the proof of Lemma 3.3) and H-B+ 1 !! = (minfgf |a(t)|) _1 = (c — 2)" 1 ~ 0.9207 (note that 
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B + is upper-triangular, whence its inverse is the Toeplitz operator with symbol a(i) _1 ), we get 
that limsup n H^" 1 )! is bounded above by (s — 2) _1 s» 2.8539. In general, it takes very large 
random matrices A n to see 1 1 ^4_~ 1 1 1 come close to sup n H-A^H because this requires a particular 
(and usually long) pattern somewhere on the diagonal(s) of A n . The latter is reminiscent of the 
finite but very long time it takes a monkey to type the complete works of Shakespeare [52]. 
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